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NUMERICAL SOLUTION OF DIFFERENTIAL EQUATIONS 
USING OBRECHKOFF CORRECTOR FORMULAS 

SUMMARY 


In this report we are concerned with the solution of differential 
equations by predictor-corrector formulas which give high-order accuracy. 

The corrector formulas are obtained from an extension of Obrechkoff' s formula. 
These formulas are characterized by the fact that they contain higher derivatives 
of the unknown function y(x) , and if their derivatives are readily determinable, 
fewer terms are needed for a given order as compared to other multistep 
methods which do not involve derivatives. Corresponding to these corrector 
formulas, certain new predictor, or extrapolation, formulas are suggested. 
These formulas also make use of the derivatives used in the corrector formulas. 


INTRODUCTION 


A class of predictor -corrector formulas of higher-order accuracy for 
the initial value problem y 1 = F(x, y) is presented here. The corrector 
formulas are based on results in a recent article [ 1] which gives an extension 
of a formula associated with the name of Obrechkoff [2] . The Obrechkoff 
formulas are characterized by the fact that they involve the higher derivatives 
of the unknown function y(x) . Corrector formulas with this property are not 
peculiar to the Obrechkoff formulas alone; less advantageously, corrector 
formulas can be obtained, for instance, by using the Euler-Maclaurin sum 
formula which also may be expressed in a form that involves the higher 
derivatives [2] . 

The great advantage in using Obrechkoff formulas is that generally 
correspondingly fewer terms are needed for a given order as compared to 
other formulas used in multistep methods which do not involve derivatives. 
This advantage is beneficial only if it is relatively easy to determine some of 
the higher-order derivatives. In some classes of problems this may be done 
by obtaining recursive formulas for the derivatives. For instance, this is 


1. This work was done at the Computer Sciences Corporation under Contract 
NAS8-18405 for the NASA Computer Center. 



the basis of the one-step method that Fehlberg [3] developed in which he 
combined a power series expansion for the solution with a Runge-Kutta process. 

We should also remember another advantage characteristic of 
multistep predictor-corrector methods; that is, the difference between 
prediction and correction gives a reasonable estimate of the local truncation 
error which can be used for automatic step-size control. 


THE EXTENDED FORMULA OF OBRECHKOFF 


Let r, m, n be nonnegative integers, r^m, Letf(x) e C m+n [a,b] 

and suppose f^ m+nf ^(x) exists and is R-integrable in [a,b] , then for every 
x in [a,b] we have 
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-f(a) = Yi -I” <x-a)V J, (a) 
j=i J * 
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so that x - a = h, then if 


If we set y = f (x) , a = x , and x = x_^ 

J n rtf-i 

f(x ) , equation ( 1) becomes 
s 


n+ 


m J M v r+ n , iJ 
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y nt^l + R ( m > n ’ r ) 
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where A., B., and R are given by equations (2) , (3) , and (4) . 


For r = 0, equation (5) reduces to the formula of Obrechkoff, thus: 
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Whereas Taylor's expansion gives f(x n ^) i n terms of f^ (x^) , 

v = 0, 1 m , the above formula of Obrechkoff gives f(x ) i n terms 

of f ^ (x^) , 0^ v — m, and / ^ ^ x n +i^ » t — v — n. Thus, the formula 

of Obrechkoff may be looked upon as a two -point Taylor type expression of 
f(x [4,5] using the values of the function and its m derivatives at x^ 

and its first n derivatives at x , . . As a matter of fact, when n = 0, 

n+i 

Obrechkoff's formula, equation (6) , reduces to the usual Taylor expansion. 
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The special case where m = n in equation (6) deserves consideration. 
From equations ( 7) and ( 8) , we get 


* _ T5 _ m - (2m - j) [ 

j j 2 m'. (m-j)l 


( 10 ) 


so that by equation (6) 
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whereby equation (9) 
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( 12 ) 


For m = 2 and 3, as examples, equation ( 11) gives 
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With appropriate predictor formulas these may be used as corrector 
formulas. For the application of equations ( 13) and ( 14) as predictor formulas 
in an n-body program see Reference 6. 


PREDICTOR FORMULAS 


Corresponding to each corrector formula we must find an appropriate 
predictor, or extrapolation, formula. Thus, as an example, in the corrector 
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* t T 

formula ( 13) for y^^, the unknown quantities y g+ ^ and y^^ appear on 

the right and must be estimated from an appropriate extrapolation formula. 
Thus, in connection with this corrector formula ( 13) , we can use 


s+i 


= y s-i + 


- 2h 2 . .. 

2h Vl + T (2y s + 


n . , 2h 5 V ... 

W + IT y (l) 


( 15) 


which we shall derive in order to illustrate how we can find predictor formulas 

to correspond to any Obrechkoff corrector formula. With the approximate value 

of y , from equation ( 15) , we form y'(x ) = F(x . , y ) and then 
S+- X S"f“ 1 S+ X S' X 

use y’ T = dF/dx = F + F F , etc. We can, of course, use recurrence 
x y 

formula where these are feasible [3] . For closer approximations, 
successive iterations between equations ( 14) and ( 15) may be made. 

To obtain equation (15) we will use the method of undetermined 
coefficients. We anticipate the form of the desired formula using undetermined 
coefficients and then impose appropriate conditions on the formula which lead 
to a system of linear equations in the solvable coefficients. Thus, we start 
with 


i = J y'(x) dx = h(a 0 + oq y'^) 


s-h 


(16) 


+ h 2 (/3 0 y^ + fa y g _ 1 ) + R , 


and it is our purpose to obtain the coefficients and remainder in equation ( 15) . 

We suppose that the remainder R is zero when y = x 11 , n = 0, 1, 2, 3, and 4 
(hence, when y(x) is a polynomial of degree — 4) . A simplification without 
a loss in generality follows if we take x^ = 0 and h = 1. We get 

a 0 + aq = 2 

- (/3 0 + dl) = 0 ( 17) 

3aq - 6di = 2 
(*1 - 3/3i = 0 
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An interpolation formula [7] is of interest here because, like the 
corrector formulas of Obrechkoff, it will give predictor formulas not only 
involving the functional values but as many derivatives as may be desired. 
When no derivatives are used, the formula reduces to the well-known Lagrange 
formula. The formulas are given by the following: 

Let f(x), together with its first m+n+1 derivatives, be continuous in 
an interval containing x 0 , x ls ... ,x n> where x.* x^. We define 


r (m + ji - i) '.ml 
c. = “ 


'i (m + n) i(m - i) !i'. 
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S k = f^) + X c! f' ^U k ) (x-Xj^.) 1 , k= 0,1,2, 
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. m+nf 1 

R - ( l) n m,n> X - f(iu+n+l) . 

R k ( ( m + n) l(m + 1) 1 f ( V 


where 0, lies somewhere between x and x. . 
k k 


(18) 

(19) 
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x. - x 
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J 1 
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Under the conditions stated 


( 20 ) 


f(x) = 


n 



x. - x 
J 


x. - x. 
J i 


S. + R 

l 


( 21 ) 


The interpolation formula is obtained by using the right member of equation (21) 
with R omitted, which is the error due to using the formula. 

It should be observed that the S k as defined by equation ( 18) need not 

all contain the same number of derivatives. Thus, S 0 could contain m 0 
derivatives, S t include m* derivatives, etc. However, this generalization 
does not affect the form of equation (21) . 
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Finally we wish to explain that in actual practice, it is usually possible 
to estimate the accuracy of the formula by computing one or two of the R. as 

indicated by the coefficients in equation (20) . 

If we assume equal intervals h between x^ and also 

that x = x^ , the above formulas are considerably simplified and are useful 
to our purposes. 


A few examples will be given as illustrations. For instance, a fifth- 
order, three-point formula is obtained if we take x = x 3 and n = m = 2. 
Then 


h <3y n-2’ 6y n-l +3y ’n ) 

+ T h!(3y n-2- 4y n-l +y n ) + R 


( 22 ) 


where we have set y = f(x) and 
R = Rq - 3Rj + 3 R 2 


and 


<*« •*k ,,y(v)( V ’ k= 0,1,2 

when 6, lies between Xq and x, 
k k 

For use with Obrechkoff corrector formulas we would probably prefer 
a two-point predictor formula (though this is not necessary) in which case 
we take n = 1. Thus, for x = x 2 and n = m = i, we get the third-order 
interpolation formula 

y ntl " ( ‘ y n-l + 2 V + h ( “ y n-l + 2 y D + E (23) 


where 


R = - R 0 + 2R t 
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and 


- ( x 2 -x )3 

R k =— 1ST- y ( V’ 


k= 0, 1 


where 0 1 lies between x 2 and x^. 

For x = x 3 and n = i,m = 3, we get the fifth -order corrector formula 

y - = ( -y . + 2y ) + 4" h ( -y , + 2y' ) + 4“ h 2 ( -y" 4 + 2y") 

J n+1 °n-l J n 2 J n-1 J n 2 J n-i J n 


i , Q , Ml ^ TM v „ 

+ — h 3 ( -y_ 4 + 2y n ) + R 
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(24) 


where 


R = - Rg + 2Rj 
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480 


k = 0, 1 


where 0^ lies between x 2 and x^, so that 

R - » 5 [l? 




Xq < 6>o < x 2 
Xi< 6 t < x 2 


(25) 


RECURRENCE FORMULAS FOR DERIVATIVES 


As an example of the use of recurrence formulas for the higher 
derivatives of our differential equation, we will give an illustration using the 
Kepler problem: 
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y= (-Tr)y 


(27) 


where 


r 2 = x 2 + y 2 

Make the substitution: 
k 

u = - — g- 
r 

And we obtain the following set of equations: 

X = ux 

y = u y 

r 2 = x 2 + y 2 

ur + 3ru - 0 . 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


where equation (33) is obtained by differentiating equation (29) with respect 
to t. 


We now introduce the following power series expansions: 
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y= Z Y n (t-t 0 ) n 

n=0 


(35) 
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(36) 


r = Z R (t -t 0 ) n 

n=0 


(37) 
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, etc. 


where the coefficients are X 

n 


J_ dx^ 

nl dt n 


Introducing equations (34) and (36) into equation (30) and equating 
like powers of (t - t 0 ) on both sides, we get a recurrence formula for the 
coefficients of equation (34) . Proceeding in this way, we get the following 
set of recurrence formulas 
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X U ! X , 
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A similar, but somewhat more complex, example may be found in Reference 3. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, September 3, 1969 
129-04-21-00-62 


10 



REFERENCES 


1. Meir, A., and Sharma, A.: An Extension of Ob reshkov's Formula. 

Siam Journal of Numerical Analysis, Vol. 5, No. 3, September, 1968. 

2. Hildebrand, F. B.: Introduction to Numerical Analysis. McGraw-Hill, 
New York, 1956, pp. 231 and 232. 

3. Fehlberg, E.: New One-Step Integration Methods of High -Order Accuracy 
Applied to Some Problems in Celestial Mechanics. NASA TR-248, 
October, 1966. 

4. Collatz, L. : The Numerical Treatment of Differential Equations. 
Springer-Verlag, Berlin, Third Edition, 1960, p. 19. 

5. Hummel, P. M., and Seebeck, C. L.: A Generalization of Taylor's 
Formula, The American Mathematical Monthly, Vol. LVI, No. 4, 

April, 1949. 

6. The N -Body Space Program. General Dynamics, Convair, San Diego. 

7. Hummel, P. M., and Seebeck, C. L.: A New Interpolation Formula. 

The American Mathematical Monthly, Vol. LVin, No. 6, June-July, 

1951. 


NASA-l angley, 1969 19 M527 


11 



National Aeronautics and Space Administration 
Washington, D. C. 20546 


OFFICIAL BUSINESS 


FIRST CLASS MAIL 



POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


l < 


l 4 



>0 v 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


" The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
■ of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



